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Elliptic flow measurements at RHIC suggest that quark gluon plasma flows with very little vis- 
cosity compared to weak coupling expectations, challenging theorists to explain why this fluid is so 
nearly "perfect" . It is therefore vital to flnd quantitative experimental information on the viscosity 
of the plasma. We propose that measurements of transverse momentum fluctuations can be used to 
determine the shear viscosity. We use current data to estimate the viscosity-to-entropy ratio in the 
range from 0.08 to 0.3, and discuss how future measurements can reduce this uncertainty. 
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Measurements of elliptic and radial flow at RHIC 
are described by viscosity-free hydrodynamics, indicating 
that the quark-gluon system produced in these collisions 
is a nearly perfect liquid QQ, 0, In particular, the 
strong suppression of flow due to shear viscosity predicted 
by weak-coupling transport calculations is not observed 
This result is exciting because a small viscosity rel- 
ative to the entropy density of the system may indicate 
that the system is more strongly coupled than expected: 
The collisional shear viscosity is proportional to the mean 
free path, which is shorter when the coupling is stronger. 
But is the viscosity really small? Hirano et al. point out 
that color glass condensate formation mayproduce more 
elliptic flow than considered in refs. P, ]M, requiring a 
larger viscosity for agreement with data j5(. 

We seek an experimental probe of viscosity that is in- 
dependent of elliptic flow. To that end, we propose that 
transverse momentum correlation measurements can be 
used to extract information on the kinematic viscosity, 

jy = rj/Ts, (1) 

where rj is the shear viscosity, s is the entropy density 
and T is the temperature. This ratio characterizes the 
strength of the viscous force relative to the fluid's inertia 
and, consequently, determines the effect of 77 on the flow 
[3 . We argue that viscous diffusion broadens the rapidity 
dependence of transverse momentum correlations, and 
then show how these correlations can be extracted from 
measurements of event-by-event pt fluctuations. 

A number of experiments have studied transverse mo- 
mentum fluctuations at SPS and RHIC 0,0- Interest- 
ingly, the STAR collaboration reports a 60% increase of 
the relative rapidity width for pt fluctuations when cen- 
trality is increased While the STAR analysis differs 
from the one we propose, model assumptions provide a 
tantalizing hint that the viscosity is small. 

Any experimental information on the kinematic viscos- 
ity of high energy density matter is vital for understand- 
ing the strongly interacting quark gluon plasma. The- 



orists had long anticipated a large collisional viscosity 
based on weak coupling QCD ^ and hadronic computa- 
tions jlCll |. with values of rj/s roughly of order unity for 
both phases near the crossover temperature ~ 170 MeV. 
Supersymmetic Yang Mills calculations give the signifi- 
cantly smaller ratio rj/s = l/An in the strong coupling 
limit . Lattice QCD calculations of the shear viscosity 
will eventually settle the question of the size of the viscos- 
ity near equilibrium 12]. However, the effective viscos- 
ity in the nonequilibrium ion-collision system may differ 
from these calculations. In particular, plasma-instability 
contributions can also explain the small viscosities in nu- 
clear collisions 13] . 

We begin by formulating a simple model to illustrate 
how shear viscosity attenuates correlations due to fluctu- 
ations of the radial flow. Next, we show how transverse 
momentum fluctuations can be used to measure these 
correlations. We then demonstrate the impact of viscos- 
ity on the rapidity distribution of fluctuations. Finally, 
we explore the implications of current fluctuation data. 

Before wading into the quark-gluon liquid, it is useful 
to recall how shear viscosity affects the flow of more com- 
mon fluids. In a classic example of shear flow, a liquid 
is trapped between two parallel plates in the xy plane, 
while one plate moves at constant speed in the x direc- 
tion. The fluid is pulled along with the plate, so that v^. 
varies with the normal distance z. The viscous contribu- 
tion to the stress energy tensor is then 

T., = -vdvjdz; (2) 

see ref. for a general treatment. 

Central nuclear collisions produce a high energy den- 
sity fluid that flows outward with an average radial ve- 
locity Vr- In the hydrodynamic description of these col- 
lisions, we typically assume that varies smoothly with 
spacetime (t, x) and is the same for all collisions of a fixed 
impact parameter. For central collisions, Vr is radially 
symmetric. More realistically, small deviations u(x) of 
the radial flow occur throughout the fluid, varying with 
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each ion-collision event. Such deviations occur, e.g., be- 
cause the number and location of nucleon-nucleon sub- 
collisions varies in each event. 

Viscous friction arises as neighboring fluid elements 
flow past each other. This friction reduces u, driving 
the velocity toward the local average Vr- The final size 
of the velocity increment u depends on the magnitude of 
the viscosity and the lifetime of the fluid. 

In order to illustrate how the damping of radial flow 
fluctuations depends on the viscosity of the fluid, we in- 
troduce a velocity increment in the radial direction u 
that depends only on the longitudinal coordinate z and 
t. Our aim is to determine the linear response of the fluid 
to this perturbation. For simplicity, we take the unper- 
turbed flow as slowly varying, and work in a co-moving 
frame where Vr locally vanishes. As in the flow of 
neighboring fluid elements at different radial speeds u{z) 
produces a shear stress 



-rfdu/dz. 



(3) 



This stress changes the radial momentum current of the 
fluid, which is generally Tq^ — 7'^(e + p)vr for energy 
density e, pressure p, and 7 = (1 — The 
perturbation u results in the change (7t(x) = STor « (e -I- 
p)u in the co-moving frame. On the other hand, the 
energy-momentum conservation law d^jT^''' — implies 
dgt/dt = -dT,r/dz. 

We combine these results to obtain a diffusion equation 
for the momentum current 
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to linear order, where the kinematic viscosity is given 
by , since e + p ^ Ts for small net baryon density. 
Observe that (^J applies for any fluctuation gf for which 
V • gt = our physically-motivated radial gt{z,t) is 
a specific instance of such a flow. Such shear modes are 
related to sound waves (compression modes) but diffuse 
rather than propagate. Note that the scale over which 
sound is attenuated Tg = (4ri/3 + C)/Ts depends on both 
shear and bulk viscosity [ij, ll5j . 

Viscosity tends to reduce fluctuations by distributing 
the excess momentum density gt over the collision vol- 
ume. This effect broadens the rapidity profllc of fluc- 
tuations. We write Q in terms of the spatial rapidity 
y = 1/2 \n{t + z) / (t — z) and proper time r = (t^ -z^)^/^ 
to flnd dgt/dr = {v / T^)d^ gt / dy^ . A similar equation is 
used to study net charge diffusion in ref . jl6l | , and we can 
translate many of those results to the present context. 
Deflning V = {{y - {y)f) = J y'^gtdy/ J gtdy for (y) = 0, 
we compute the rapidity broadening 
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FIG. 1: Rapidity spread vs. time for momentum diffusion 
computed using (|^ and ||7|l for the large viscosity (viscous) 
and small viscosity (perfect) scenarios discussed in the text. 
The gray area marks the range extrapolated from data in 
ref. 1^ using lHH . 



We extend this discussion to address a more general 
ensemble of fluctuations by considering the correlation 
function 



(3t(xi)5t(x2)) - (gt(xi))(.gf(x2)). 



(6) 



g^eq- The spatial 
eq is broadened by 



In local equilibrium, Vg has the value 
rapidity dependence of Avg = 
momentum diffusion. If the rapidity width of the one- 
body density follows jSJ, then the width of Ar^ in the 
relative rapidity j/r = J/i ^ 2/2 grows from an initial value 
CTo following 
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where Tf is the proper time at which freeze out occurs. 
This equation is entirely plausible, since diffusion spreads 
the rapidity of each particle in a given pair with a vari- 
ance AV. We then take 



Ar,(2/,,,y„)oce-^^/2-^-^»/2^^ 



(8) 



where AV = V — V{to) for Tq the formation time. 



where Q gives the width in relative rapidity and the 
width in avera ge r apidity ya — {yi + y2)/2 is S. We 
assume S » cr [igl- Observe that (0) and ((Sj) are exact 
for our diffusion model 

Gyulassy and Hirano survey computations of the ratio 
of the shear viscosity to the entropy and flnd that both 
the hadron gas and the perturbative quark gluon plasma 
have rj/s 1, if one naively extrapolates these calcula- 
tions near Tc These values correspond to — rj/Ts 
roughly of order 1 fm for Tc = 170 MeV. On the other 
hand, they argue that the entropy increase near Tc re- 
duces rj/s for a strongly interacting plasma, perhaps to 
the supersymmetric Yang- Mills value rj/s ~ I/Att. 

Motivated by these estimates, we show the increase of 
a given by © and as a function of proper time r for 
two extreme and highly schematic scenarios in flg. H In 
the 'perfect' scenario, we take v ~ 0.1 fm for the plasma 
and mixed phase, and v ^ 1 fm for the hadronic phase. 
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In the 'viscous' scenario, we take i/ ~ 1 fm for the entire 
evolution. In both cases, we assume that the formation, 
hadronization, and freeze out times are 1 fm, 9 fm, and 
20 fm respectively. 

We stress that the rapidity width depends on the vis- 
cous diffusion coefficient integrated over the collision life- 
time. Comparing the viscous and perfect scenarios in 
fig- Q we see that the largest contribution to this width 
comes from the earliest times. Consequently, we expect 
measurements of this width to yield information on the 
viscosity when the evolution is dominated by partons. 

Variation of the radial fluid velocity over the collision 
volume induces correlations in the transverse momenta 
Pt of particles To describe such correlations, we di- 
vide the inhomogeneous fluid into cells small enough to 
be uniform. Particles emerging from cells of different ra- 
dial velocity Vr are more likely to have different pt than 
particles from the same cell. The number of particles 
of momentum p in a cell at position x at the instant of 
freeze out is dn = f{x.,p)dpdx, where dp = d^p/{2-K)^ 
and dx = d^x. We take /(x, p) to be a Boltzmann dis- 
tribution corresponding to a fluid velocity v(x) and a 
temperature T(x) that vary with each event. A similar 
formulation is used in ref. 18] to compute nonequilib- 
rium Pt fluctuations. Here, we focus on central collisions 
where local equilibrium is likely achieved. 

To characterize the dynamic correlations of pt , we use 
the transverse momentum covariance 

C^{N)-^Y.PuPt,)-(pt)\ (9) 

where i labels particles from each event and the brackets 
represent the event average. The average transverse mo- 
mentum is {ptj = (^Pti)/{N). This covariance vanishes 
in local equilibrium, where the momenta are uncorrelated 
and number fluctuations satisfy Poisson statistics. 

This covariance is related to the spatial correlations of 
the momentum current © by 

C = (N)-^ J Arg{xi,X2)dxidx2. (10) 

To obtain this result, observe that {N){pt) = (/ ptdn^ = 
J {gt{:x.))dx, where 

5tW = J fiyi,p)ptdp (11) 

is the momentum current discussed earlier. Similarly, we 
write the unrestricted sum ^puPtj = J PtiPt2dnidn2 ~ 
J gt{pci)gt{x.2)dxidx2 and average over events to find 

J rgdxidx2 = ( PuPtj) - {N)'^{pty 

all i,j 

= {Nfc + iY^pD; (12) 
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FIG. 2: Rapidity dependence of the pt covariance @ for the 
scenarios in fig. 1. The initial distribution has ao ~ 0.5. 

the second equality follows from Q . In local equilibrium, 
C = implies J Tg cqdxidx2 = (^Pti)- Subtracting this 
term from H12|) gives (|10|l . 

The correlation information probed by C differs from 
that found in the multiplicity variance R = {{N'^) — 
(Nf - {N))/{Nf [13. As before, we write R = 
(N)^'^ J Arndxidx2, where Ar„ — — fn, oq and 

r-n = (n(xi)?i(x2)) - (n(xi))(n(x2)). (13) 

The density correlation function 113|l carries different in- 
formation than 10) because particle number is not con- 
served. Density ffuctuations evolve by the full hydrody- 
namic equations, while gt follows diffusion. The corre- 
lation function probed by net charge fluctuations is dis- 
cussed in ref. [la| ■ 

Viscosity information can be obtained from C as fol- 
lows. The broadening in rapidity of Avg depends on the 
shear viscosity via 0. Equation (|10|l implies that the ra- 
pidity dependence of Avg can be measured by studying 
the dependence of Q on the rapidity window in which 
particles are measured. We illustrate this acceptance de- 
pendence in fig. 2 for our idealized scenarios by integrat- 
ing © over the interval -A/2 < yi, y2 < A/2; {N)C^ 
is the value for large A. 

The STAR analysis in ref. ^ incorporates some of 
these ideas and, intriguingly, finds a broadening in ra- 
pidity together with a narrowing in azimuth for pt corre- 
lations in central compared to peripheral collisions. We 
will use the rapidity information to estimate the viscosity. 
However, the measured quantities differ sufficiently from 
C that this estimate requires significant model assump- 
tions. We therefore regard the result only as a signal of 
our method's promise. 

STAR employs the transverse momentum fluctuation 
observable Aa^^ to construct a correlation function as a 
function of rapidity and azimuthal angle. They find that 
near-side correlations in azimuth are broadened in rela- 
tive rapidity, with a rapidity width ct* that increases from 
roughly 0.45 in the most peripheral collisions to 0.75 in 
central ones We estimate AalJ {N) C - {pt)^R. 
This quantity therefore depends on both momentum cur- 
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rent and density correlation functions ^ and (|13|l . 
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j {Arg - (pt)2Ar„}dy] 



dv2 



(14) 



(T^ p for the 



We can directly compare cr* to a in fig. 1 if Ar^ and 
Ar„ have the same widths. Equation Q then implies 
that the widths in central and peripheral collisions satisfy 
'^c ~ '^p = 4i'(t^p — c)- Taking the freeze out times 
in a central and peripheral collisions to be Tp^ c ~ 20 fm 
and tf,p ^ 1 fm, we then find v ~ 0.09 fm. The value 
TF, p ~ 1 fni is reasonable, since ref. |^ argues that the 
average participant path length is about 1 fm for these 
peripheral collisions. We use ^ to find rj/s ^ 0.08. 

This result is remarkably close to the supersymmetric 
Yang Mills value l/47r, and is consistent with some hydro- 
dynamic comparisons to elliptic flow data 0|- However, 
we must be cautious: If Arg and Ar„ have different ra- 
pidity widths a and (t„ then their relation to cr* depends 
on the relative strength of these contributions. Data in 
ref. [lOl may indicate that an is roughly twice cr*. Gen- 
erally, a is bounded by cr„ and cr*, since ()14|l implies 
ctJ w cr^ + P[a'^ — <^n)- Although f3 is not measured, the 
largest a may be is an ~ 2(7*. For the maximum value 
a = 2(7*, our dynamic assumptions yield rj/s = 0.3. To- 
gether, our estimates constitute an uncertainty range for 
the viscosity-to-entropy ratio, 0.08 < rj/s < 0.3 

In fig. n we indicate the range of a^ ^ 
most central and peripheral STAR values from ref. 
The gray band follows from the uncertainty in relating 
(7 to (7*, i.e., (7* < (7 < 2(7*, which greatly exceeds the 
experimental uncertainty. Our 'perfect' liquid curve falls 
above the bottom of this range because we assume a large 
viscosity following hadronization. 

In summary, we find that shear viscosity can broaden 
the rapidity correlations of the momentum current. This 
broadening can be observed by measuring the trans- 
verse momentum covariance ((HJ as a function of rapid- 
ity acceptance. Our rough estimate from current data, 
77/s ~ 0.08 — 0.3, is small compared to perturbative com- 
putationsi. To reduce the uncertainty range, we sug- 
gest comparing C to allow more direct access to the mo- 
mentum density correlation function. In principle, freeze 
out times can be inferred from other measurements |2ll |. 
Note that minijets, color glass, and other contributions 
to the particle production mechanism influence the initial 
fluctuation spectrum and, correspondingly, modify ao in 
(0 . We assume that this contribution cancels in studying 
the centrality dependence at a fixed beam energy. 

The viscosity of a common fluid can be measured by 
applying a known pressure and observing the resulting 
flow in a fixed geometry, e.g., a pipe. Alternatively, one 
can study the attenuation of high frequency sound waves 
from a calibrated source. Efforts to compare radial and 
elliptic flow measurements to viscous hydrocivnamic cal- 
culations are analogous to the first method [31 ■ Our ob- 
servable C is in the spirit of ultrasonic attenuation. The 



early dynamics produces a spectrum of fluctuations anal- 
ogous to sound waves that are attenuated by viscosity. 
We suggest that experimenters pursue both approaches 
to extract quantitative viscosity information from ion col- 
lisions, since the geometry, initial conditions, and probe 
parameters are all unknown. 
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